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Abstract 

We determine a considerable class of nonlinear partial differential equa¬ 
tion systems, related to a large class of fluid mechanical models, which 
have global regular solutions. Uniqueness is not a direct general conse¬ 
quence of this method. The solution scheme uses local-time represen¬ 
tations in terms of probability densities and applies in strong function 
spaces. We observe that strong Gaussian damping of scaled models can 
offset possible growth of subhomogeneous nonlinear terms, where the spa¬ 
tial scaling parameter r > 1 indicates the deviation from a strong semi¬ 
group property in strong norms. The application to the Navier Stokes 
equation is considered as an example of the scheme in detail. Although 
this construction of global solution branches works for a large class of 
models (and can be extended to models with variable viscosity terms and 
even to models with highly degenerate diffusions), we remark that a direct 
application to randomized models (for example average models or models 
with white noise) is not possible in general, as energy can be transported 
to higher frequencies caused by stochastic mixing. We also compare with 
results obtained from singularity analysis and CKN theory, and consider 
generalisations to models with variable diffusions or highly degenerate 
diffusions. 
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1 Definition of a integrable class of nonlinear 
partial differential equation systems 

For a positive viscosity constant zz > 0 consider the Cauchy problem 

Vi^t - v^Vi + fi{v, Vv) = 0, 1 < i < D, (1) 

where v = (vi, ■ ■ ■ , vb)^ and Vv = (Vui, • • • , Vvd)^ are defined on the domain 
[0,oo) X with data 

?;i(0 ,.) = hi, 1 <i < D. (2) 

Here, D > 3 denotes the dimension of the problem and R. denotes the field of 
real numbers. 
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Remark 1.1. The operatorsin m can be global operators, e.g. partial integro- 
differential equation as in the case of the(Leray projection form of the) incom¬ 
pressible Navier Stokes equation. We shall define a scheme where a solution is 
bounded in a ball of finite radius C > 0 with respect to a regular norm, and 
such that this bound is preserved in time, i.e., there is a finite constant C > 0 
such that for m > 2, all tg > 0 and some A > 0 (independent of to > 0) we have 

maxi<i<£) \ < C, and 


maxi<,;<_D Olz/mnc™ ^ maxi<i<_D \vi{to + < C. 

(3) 

In general our construction implies that C > maxi<,;<£) \ and a spa¬ 

tial scaling parameter indicates this deviation of a strong semigroup property. 

Assume that the following conditions are satisfied. 

a) For m > 2 and 1 <i < D the data satisfy 

h, (4) 

Here, C™ is the function space of continuous functions with continuous 
multivariate derivatives up to order m, and is the standard Sobolev 
space of order m, i.e., the space of functions with multivariate weak deriva¬ 
tives in up to order m. For m = 0 the function space C := (7° is the 
set of continuous functions. We shall use the norm 

0<|Q!|<m 

where for a continuous and bounded function g : K.^ —> K. we define 

kL2nc ■= I5L2, sup |g(a;)| \ . ( 6 ) 

In some situations, it is possible to transform to problems on compact 
domains. If H C K." is compact, then |5|^2(f2)n(7(f2) denotes a local form 
of this norm. 


b) The nonlinear terms fi satisfy a sub-homogeneity condition in the sense 
that for spatial scaling y = rx, r > 0 (where r is any given positive real 
number), and for functions Vi, 1 < i < D with Vi{t ,.) £ m > 2, 

we have for all 1 < i < T>, all t > 0, and all y = rx € 


\fi{v{t,x),Wv{t,x))\ < r fi{v{t,x),-Wv{t,x)) 

I I j. 


(7) 


Note that for the transformation vl{t,y) = Vi{t,x) we have 


rf, ( v{t,x),-Vv{t,x) ) = rfi{v^{t,y),Vv^{t,y)), 


( 8 ) 


because Vij = Hence, the inequality in ([HI) may be also 

rewritten such that for any given r > 0 and with vl{t,y) = Vi(t,x) we 
have 


\f^{v{t,x),Wv{t,x))\ < r fi{v^{t,y),Wv''{t,y)) 


(9) 
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c) In the following for functions F'^ we use the notation 

•= Vu(t,a:)), 

i.e., the notation refers to derivatives with respect to the argument Xj of 
the composition; some readers may prefer the notation 

■^F^{v{t,x),Vv{t,x)) = {F^iv{t,x),Vv(t,x))) 

which we avoid for simplicity of notation. Multivariate derivatives with 
respect to the arguments of Fj may be denoted by where in 

this case 7 is a inultiindex of length D + D^. Note that these are usual 
multivariate derivatives. We assume that there exists a matrix of regular 
functions {Fj)i<ij<D such that for all 1 < i < D 

D 

^F;,^-KV7;) = /,(u,Vr;). ( 10 ) 

i=i 

Remark 1.2. The assumption of a matrix F^ with m implies that we have 
classical local time representations in terms of the first spatial derivatives 
of the Gaussian. Symmetries of the first order spatial derivatives of the 
Gaussian can be exploited if the convoluted term in the local represen¬ 
tation satisfies a Lipschitz condition. This structure together with the 
scaling assumption of item b) ensure that local diffusion damping can 
offset possible growth of the nonlinear terms for a spatially scaled model. 

We assume that fi satisfies a local Lipschitz condition, i.e., for g = 
{91 r " }9 d), 9i G G"* n iL'” and for any finite constant G > 0 there 
exists a finite constant L (depending on G) such that 

Note the loss of one order of regularity in this stipulation. Further¬ 
more, if Fj has a local interpretation then we assume that for all mul¬ 
tiindices 0 < |/3| < m the functions D^^^Fj : R are lo¬ 

cally Lipschitz continuous (on any compact domain in H C R^ and 
with respect to all arguments). In any case we require that for any 
9 = ( 51 , • • • ,9d) with < G < 00 we have Lipschitz continuity 

of y —>■ DJFj {g{y),yg{y)) with hnite Lipschitz constants LF (dependent 
on G) in the sense that for all 1 < i,j < D and all 0 < I 7 I < m — 1 and 
y, y' € R^ 

|Gj( 5 (y),Vy(y))-Fj(y(y'),V 5 (y'))| <Gf|y-y'|, (12) 

and 

\DlFl^{g{y),Vg{y))-DlFl^{g{y'),Vg{y'))\<L^^\y-y'\. (13) 

Here, denotes the multivariate spatial derivative with respect to the 
multiindex 7 = ( 71 , • • • , jd) and with respect to the argument of g. 


3 


d) The verification of the technical condition in c) can be simplified for spe¬ 
cific models if we add a stronger assumption concerning the data. This 
additional assumption is also useful, if generalisations of the diffusion term 
or viscosity limits are considered, e.g., if we replace the Laplacian term 
by a Hormander vector-field condition for highly degenerate operators of 
second order. Let 

/i, G if'" n c™ n (i4) 

where for a given integer I > 1 

eUm = {/ : ^ R : 3c> 0 V|x| > 1 VO < I 7 I < m \D2f{x)\ < 

(fs) 

Note that the latter function space has a multiplicative property, i.e., 
g,h G Cpgj ^ implies that gh G In addition to (fT^ we then require 

a ’submultiplicative property of order G {0, • • • , m{D -b 1) — 1}’ , i.e., 
for g = {gi,--- ,gD) along with gi G we require that for some 

k G {0, • • • , m[D -b 1) — 1} we have 


g, G e 


m(D+l) 
pol,m 1 


(VgOj e e 


m(D+l) 

pol,m 


implies i 4 fF;(g,V 5 )Ge 


2m(D+l) —fc 
pol,m 


(16) 

Even the latter condition is a sufficient for a class of operators which 
includes some models of fluid mechanics. Especially, this ’multiplicative 
property’ of the nonlinear term holds for the incompressible Navier-Stokes 
equation operator. The number m{D + 1) in the upperscript of 
is not a sharp choice, but a choice which may be also sufficient for the 
consideration of viscosity limits. 


We have 


Theorem 1.3. If the set of conditions a), b), c) or the stronger set of conditions 
a), b), c) and d) are satisfied, then the Cauchy problem in (Q)), 0) has a global 
classical solution Vi G ([0, 00 ), fl C""), 1 < * < £>. 

Some remarks are in order. 

Remark 1.4. We have described a quadratic system where a vector v with D 
components solves D equations. This is not an essential restriction, and it is 
convenient. 

Remark 1.5. We do not claim uniqueness in (IE31) although in many special 
situations standard arguments may lead to uniqueness. Note that in some situ¬ 
ations the method described below may be applied in order to get global solution 
branches of viscosity limits of the equations in o, and then there are examples, 
where determinism is lost and global regular solution branches exist next to 
singular solutions. 

Remark 1.6. We use Gaussian upper bound estimate. We cannot prove a strong 
contractive semi-group property for a class of nonlinear operators which includes 
the Navier Stokes equation operator. The deviation from a strong semigroup 
property is indicated directly, if we estimate for v^, 1 < * < 4? with r > 1 . 
In any case estimates are obtained for multiples of a discrete time, and then 
the upper bound for all time is constructed by a local time contraction result 
in regular space, which gives an upper bound in regular space, and this upper 
bound is larger than the data in the same regular space in general. 
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Remark 1.7. We may reduce to a local Lipscliitz continuity assumption in c) if we 
add d) or the weaker assumption that for u = (ui, • • • ,ud) with m G fl C™ 
for m > 2 we have that Fj{u, Vtt) € L^. 

Example 1.8. In case of the incompressible Navier Stokes equation (cf. [1] for 
the modeling) the incompressibility condition implies 


^ d{viVj) 
^ dxj 


D 


E 



D 

+^*E 

i=i 


dvj 

dxj 


D 


E 



Therefore we may define 


Fj — VjVi SijKjj ^sp 


D 

l,m—l 


(17) 


(18) 


where Ko denotes the Laplacian kernel of dimension D and 6ij is the Kronecker- 
S. Here *sp denotes the convolution with respect to the spatial variables. Note 
that for y = rx, vl(t, y) = Vi(t, x), and z = rw, and the first spatial derivative of 
the Laplacian kernel x —>■ KD^i{x) = in the Leray projection term (written 
as an operator on the Jacobian J{v) = rn<D) satisfies 

HJ{v)) = SE=1 Vl,m{t,w)Vm,l{t,w)dw 


= /rd EE=i z):^dz 


(19) 




This linear spatial scaling of the Leray projection term can be observed also 
from the linear scaling of the pressure gradient. The Poisson equation for the 
scaled equation becomes i — and is identical to the 

pressure elimination equation in original coordinates as the factor cancels. 
More explicitly, for y = rx and p^(y) = p(x) we have the spatial derivative 
transformation (V)ip := p^i = and the original Navier 

Stokes equation transforms to 


dt 


D 


^ - r^izAvl + rJ2 


1=1 


dv^ 


( 20 ) 


and, applying the divergence operator and incompressibility, i.e., divv = 0, 
leads indeed to 

D D 

E ^ E <^<,1 = ( 21 ) 

l,m—l l,m—l 

Hence, the Leray projection form of the scaled equation becomes 

o,,r p)„,r p 

^+ =r / KD,i[. - y) ^ v\ v'!;^i{y)dy (22) 
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If the additional assumption d) holds, then it is rather obvious that the technical 
Lipschitz condition in c) holds, but it can also be verified under the weaker set 
of assumptions a), b) and c). Assume that at time to > 0 we have the upper 
bound 

max IvUtn, .)| < C for some m> 2, (23) 

i<i<D I ‘ ^ no 

for some finite constant C > 0. Local time iteration schemes with respect strong 
norms lead to local time representations 


< = <(to,.) G: + r * Gl 

+r KdA- - y) <m(-> V^V * 


( 24 ) 


on a short time interval [fo 7 fo + A] for small A > 0, and where * denotes 
convolution with respect to space and time on this time interval. Here, G)) is 
the fundamental solution of q^t — = 0 (considered on the time interval 

[to,io + ^])- Similar local time representations hold for spatial derivatives of 
velocity function. Indeed for multiindices /3 = (/3i, • • • , /3d) and 7 = ( 71 , • • • , 7 ^) 
and 1 < |/3| = 1 + I 7 I < m with = 7 ^ + 1 and /3j = jj for p ^ j we have on 
the time interval [to, to + A] 

DPv: = .) GD + rD2 (Ef=i * GO,^ 

( 25 ) 

(Im^KdAv) - y)dy) *Glp. 


Local time contraction shows that 


max sup |t',’^(t,.)! < G+1 


such that we have the Leray projection estimate 


{{Kd,. *sp D2 (vAva)) * GO,,) (t,.) 


L^nc 


( 26 ) 


< rc^{G + 1) max^^ max|5|<^_i 


{{KD,^ Div\)*Gl^^) (t,.) 


L^nc’ 


(27) 

where is a finite constant depending only on dimension D and the regularity 
number m. Since for t € [to, to + A] and 0 < |i5| < to — 1 

\Div[(t,.)\^,^^,<C+l ( 28 ) 


we have spatial Lipschitz continuity of KD,i*spG>Ai {t ,.) for multivariate spatial 
derivatives of order |(5| < to — 1 and in case of dimension D = 3. As in the main 
theorem we then may use symmetry of the first order derivative of the Gaussian 
(cf. below). Indeed, application of Theorem [T|3] implies the existence of a global 
regular solution branch. In this special case uniqueness is implied. 
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2 Proof Theorem 11.3 


i) We rewrite the equation using the integrability of the nonlinear terms. 
From 0 and assumption c) we have 

D 

Vi^t — i^Avi + ^2 1 < i < D. (29) 

Recall the notation Fjj{v,\7v) = [Fj{v,Vv)) where the derivative is 
with respect to the jth spatial argument Xj of the composition of F"^ with 
{v, Vv). 

ii) Assume that data Vi(to,.), 1 < i < F are given. We have local time 
contraction in spatial iJ™ fl C™ space. More precisely, define a time local 
iteration scheme wf, 1 < i < F, fc > 0 in a time interval [to, io + A], where 
to > 0 and where for t G [to, to + A] we define 

z;°(to,.) = v^{to ,.) *sp Gl, (n.b. .) = h,{.)). (30) 

Furthermore, for k > 1 v!^ is a, solution of the linearized equation 

D 

^ FF(v'^-\Wv>^-^) =0,l<t<F, (31) 

where uf (tg,.) = Wi(to,.) for 1 < t < F. Note that 

D 

uf = ^.(to,.) Gl + Y, * Gl (32) 

1=1 


for 1 < 1 < H. 


Lemma 2.1. Let to > 0 and assume that for some m > 2 and a finite 
eonstant Co > 0 we have 


Then there exists a A > 0 dependent only on dimension and on the con¬ 
stant |fi(to, such that on the time interval [to, to + A] the func¬ 
tional increments = v^'^^ ~ ’^2 ^ — J — ^ satisfy for k > 1 


sup 

t6[io,to+A] 




H^nc^ — 


1 

< - 
2 


sup 

i6[to.io+A] 


\6v’l{t,.)\ 




(34) 


and 


sup \Svl{t,.)\^^^^,, 



(35) 


The proof uses classical representations and the Lipschitz continuity of 
(spatial derivatives of) C® and is given in the appendix. 

iii) The local time contraction result of item ii) implies that there exist regular 
local time solutions Vi £ C^ ([to, to + A], iJ"® fl C"®), I < i < F. 
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Remark 2.2. If the fi satisfy the condition in d), then for the functional 
increments 5v^ = , /c > 1, we have 

5v^ e Cl ([to, to + A], i?™ n c™ n , i < * < c, 

and this holds also in the limit for Wj, 1 < t < It. Here, note that 


= E *Gl-Y^ Fl^{v^-\Vv^-^) * Gl (36) 

i=i i=i 

for 1 < i < D, where the submuliplicative property of fi implies a preser¬ 
vation of the order of the order of spatial polynomial decay (as k increases) 
where it clearly offsets possible effects of decrease the order of spatial poly¬ 
nomial decay caused by convolutions with a first order spatial derivative of 
a Gaussian. For more complicated second order diffusions we have a con¬ 
volution with a density which can cause a stronger decrease of polynomial 
order decrease at spatial infinity, and we need a stronger submultiplicative 
property (smaller k in assumption d)). 


We consider local classical representations of solutions and use the convo¬ 
lution rule in order to get for all t £ [toj^o + A] and x G 


Vi(t,x) = Vi(t0, ■) *sp Gy + PjJ * Gy 

= Vi(to, .) *sp Gy + J2f=i P'j '^1') * Gyj 


(37) 


Again remember that , j refers to derivative with respect to Xj , where we 
use the notation {v,Vv) = (Fj^ {v,Vv)) ^ Furthermore, the symbol 
* denotes convolution with respect to space and time. In the last line the 
nonlinear terms are convoluted with first order spatial derivatives of the 
Gaussian, while the first term on the right side of (1571) is a convolution 
with the Gaussian, which behaves completely different. Here we may use 
only one spatial scaling parameter r > 0 and exploit spatial effects of the 
operator. We shall observe that small damping of the latter term offsets 
possible growth caused by the former. Note that for spatial multivariate 
derivatives of order \(3\ and for 0 < jyj -|- I = |,5| < to, (5k = Ik F 1, and 
/5i = Ji for I ^ k we have time local representations of the form 

DPv,it,x) = DPv,{to ,.) Gy + DZFl^ iv,Vv) * Gy,k. ( 38 ) 


iv) Using Lipschitz estimates we get an upper bound of the nonlinear terms. 
First note that first order spatial derivatives of the Gaussian Gyj have a 
symmetry which can be combined with Lipschitz continuity of the nonlin¬ 
ear function terms. We have 


Gyj {t,y) 


-Vj 1 


exp 



(39) 


Define y^' = {y{’ ,■■■ .y^ ), where = yk for j and yj’ = -y^. 
Assuming |ni(to, ■)|jjmpic'm = < oo the local contraction results tells 

us that for some A > 0 and all time t G [to, + A] we have 

|ni(t, OlcmnH™ - ^*0 
8 


( 40 ) 




Recall that functions Fj : —>• R are Lipschitz on any finite ball, 

especially on a ball of Bct +i(0) of radius Ctg + 1 around zero. It follows 
that for all i < i,j < D 

F;{v,yv){t,x - y)^-^^exp dy 


(/rd - y) - F^{v,Vv){t,x - y^'-))) x 

— —^ exp dy 


2iyt 


/47TL/t 


(41) 


< 2Lf 


fu 


/47riv't 


73 - exp - 




where we use assumption c). Analogous estimates hold for spatial deriva¬ 
tives with the related Lipschitz constants L^. 

v) The nonlinear upper bound of item iv) has a scaling which is different 
from a normal Gaussian. Moreover the linear and the nonlinear part of the 
equation have scaling constraints due to the sub-homogeneity condition. 


(42) 


For t 

= rp and 

rx = y, and 

<’’’(Ay) = 

yP’"'^ 

PJ 

r, and Vijj ■ 



<■; - 

pr^vAyP"'' 4 

-Ph{v^'\r 

Now, 





\MvP 

''"{t,y),rVvP- 

II 


< r\fi{v{t,x), iVu(t,a;))| = r|/,(u’'(t, y), Vu’'(t, y)) I. 


(43) 


Hence, 


D 


, rVu^’’’) =0, 1 < * < 

t=i 


where 


= \pMvP^'',rVvP’n\ 


< pr\MvP’^,VvP’n\ = pr\Ef=iFlji'pP’^^VvP’% 

The local contraction result transfers to the scaled situation. It follows 
that on a time interval [to, to + ^o] with Aq = pA we have 


(44) 


(45) 


sup 

[ip ,to+^o] 




H’^nc^ 


i{to, .) 


H^nC’^ 


+ 1 := C, 


tom 


1 . 


(46) 


Note that this holds for to = 0 especially. 

For the scaled function we have the local representation 

x) = <’’’(to,.) GPX + pr Fl^ (vPX, VvPX) * Qp^ 


(47) 


= <’'-(to,.) GPx + pr F; (vPx, 


9 











For the scaled Gaussian G(^’[(r,x) := G,y{t,y) we have 




C;;^ _ ^ 

— {A'Ki'pr'^r)^\y\ 


-‘2yi 




. exD ( _ 

ipr^iyr J 


2)S-D/2 / \y\^ \ 

J I 47t pr'^ur J 


D/2+1-S 


exp - 


\ 4pr^iyT J 


( 48 ) 


Hence we have for 5 € (0,1) and all p, r > 0 

I (t, y) I < —- - — ^ 

' ’ ' ['i'iTpr^vTY\y\ 


S\„i\D+l-2S ’ 


(49) 


where the upper bound constant 

C= sup (z)^/^+^”^exp(-z2) >0 
hl>o 

is sufficient and independent of > 0. We are interesting in the scaling 
of convolutions of these Gaussians with (globally) Lipschitz continuous 
functions y —t /(y) with upper bound lo\y\- Now, as z/ > 0 we may choose 
the spatial parameter r to be large enough such that 

Airpr^i/ > 1. (50) 

Using the abbreviations ar = t — a, and B = {y||y| < Anpr^v}, find 
/tI/B^o|y||G;(;[(cr^,y)|dydcr < 

ftl Ib ( 4 JLJ ' ^ exp (- 4 ^S^) dyda 

^ /to Sb Cdyda < IoC*{t - toY-^Yni^pr^Y 

(51) 

for any S G (0,1) and for some finite constant C* which is independent 
of the parameters p, r,(radial coordinates with measure r^~^dr may 
be used for the latter observation). Note that for small time and/or 
small parameters the corresponding integral on the complementary do¬ 
main ■ ■ ■ becomes small such that (ICTl) is indeed the essential 

estimate in the sense that it gives the whole upper bound up to an ’small 
e’. Now we can estimate the nonlinear term rephrasing m above and 
including time. Using the sub- homogeneity property we have the upper 


10 















bound 


z)x 


-2zi 




TT exp - 


( _^ 

\ 4:upr^{T — s) I 


dzds 


Ef=i pr (f,l f^n ^vn{s, y-z)- F^(v\ Wv^s, y 


\y.i\ 


2772/5 


-y/ 47 Zu{t — s) 


exp dyds 


<prEU^Lf JIJ, 


J Jto ■ 


y/4^pr^l{r-^ ®^P ( 4pr2i.(T-s) 


i»r 


)) ‘^y 



X 


< prEf=i 2Lf{ATTpr^v)\C*{T - to)^-^ + e), 

where we have a finite constant C* with 

C*{t - to)^-^ + e = 2D ^^j^i§+T^dyda 

+ Ef=i Itl /r^>\bi(o) \yj\\GZ(<^’y)\‘^<^dy. 


(52) 


(53) 


Here e is defined by the second term on the right side of (IMl) . The first 
summand on the right side of (1531) is independent of r, p, v. The last 
summand on the right side of (15311 is relatively small (goes to zero with 
exponential decay as pr^v or Aq go to zero). Hence e is small compared 
to any of the polynomials of lower order which determine the main part 
of C*. Summing up we have 


Ef=i pr ftl fxn FJiv’', Vu’-)(s, y - z)- 


Avs 


-y/47r2/(r —s) 


exp(- 4.pi"[!-s) ) dz 


< prLQ{AT:pr‘^vf{{T-t^Y + e), 


where Lq := (X)j=i 2T*°) and with an e > 0 of exponential decay with 
respect to i/pr^Ag as i/pr^Ap becomes small. Analogous considerations 
lead to an upper bound 


Ef=i pr ftl fxn Fdpjiv’', VvY(T-s,y-z): 


■ exp - 


f _ 

\ Ai'pr^t J 


dz 


< prlmYnpr'^izY {{t - toY ^C* + e), 

(55) 

for some constant Lm which is independent of p, r and u and which we 
choose such that it serves for all spatial derivatives up to order m. We 
may assume that Lm A 0 such that we may use Lm, for all these upper 
bounds as a constant which is independent of specific multiindices less or 
equal to m. 

vi) We use small damping of the heat convolution on infinite space, if a thresh¬ 
old is exceeded, i.e., we may assume that for the data we have for all 
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0 < |/3| < TO 


max \D^Vi(to,-) 
l<i<D I ' 


L^nc - 


(56) 


If the latter condition is not satisfied for some (3 then there is a A > 0 such 
that the respected norm is less or equal to 1 for some time t £ [to, + A], 
and we need no damping estimate for this part of the i?™ fl C^-norm 
in the interval [to, to + A]. This way we construct an upper bound close 
to a constant Cm, where Cm is maxi<i<£) plus the number 

of terms in the standard definition of the H C'^-norm. We consider 
L^-estimates. We apply a Fourier transform with respect to the spatial 
variables, i.e., the operation 


T(u)(r, ^)= / exp (— 27 rta:^) m(t, a;)(ia 

Jro 


(57) 


in order to analyze the viscosity damping encoded in the first term on the 
right side of (1^ on a time interval [to,io + Ao], where Aq = pA. For 
T € [to, ^0 + Ao] and parameters r, p > 0 we have 


T(z;f’’'’"(to,.) Gr(r - to)) = ^(^^"(to, ■)) - to, ■)) 

= 3^ {to, ■)) exp (-47r2pr2i/(r - to)(.)^) , 
where we use (let to = 0 for simplicity) 

J(cr(r,.)) (r,0 = 3- {-^^^exp {-SC )) (P.fl 

= exp (—47r^pr^z3r|^|^) . 


(58) 


(59) 


In the following we let to = 0 and remark that the following estimates 
hold for to > 0 if T is replaced by r — to. For A > 0 small enough (such 
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that, say, STr^pr^z/rAg < 1), and for r £ [0, Ag] we get 


H'"’’’'{to, ■) *sp GP’^{t,.)\ 


L2 


= Ik^ (^ 0 , •)) (f) exp (-dTT^pr^z/rlCp))^ d( 

= /rd (^ 0 , O)^ (0 exp (-STT^pr^z/rieP)) 

(j (to, O)^ (0 exp (-STT^pr^i/rl^p)^ 
(j (to, 0)^ (0 exp (-STT^pr^z/rl^p)^ d^ 


= /i 


«^\{lijl<^,i<y<-D} 


+ f< 


{kd<A, l<t<-D} 


< 


/r^> (^ 0 , 0)^ (C) exp (STT^r^lyrA^)'j d£, 


L 


{|«d<A, l<j<D} 


(j(vr''(to,.)f(ox 


X (exp (—87r^pr^izr|^p) — exp (—STr^pr^izrAg)) jd^ 


< \J(vrn(to, oil. exp (-S^Vr^z^rAg) + (8D7r^pr^:.rA^,+^} . 

(60) 

Here, we use the assumption that Aq > 0 is small enough such that the 
upper bound estimate is a straightforward consequence of a Taylor formula 
(we may choose Aq small and Sir^pr^iyrAo < 1) and use the abbreviation 


:= sup |j(i:)fuf’’'’‘'(to,0)^(O|- 
{|«i|<A} 


(61) 


which is a finite constant (since |ui‘'(to, Ol^f.piC'. is finite). 

If we take the square root we may use the asymptotics \/l + a = 1 + + 

O(a^). 

For T £ [0, Aq] and 


0 < Aq < max 


_^_,i\ 

SzT^r^z/max{c^, 1} ’ 2 j 


(62) 


we get (the generous) estimate 

\'>J?’''’''{to, ■) *sp G^{t, .)|^2 < \3'{vl’''){to, .)|^2 exp (-dTr^i/pr-^rAg) 


+c^ {SDtt^ pr^ i/tA^'^^) 


< .)|^2 exp (-dTr^i/prVA^) + {SDtt^ pr^iyTAl+^) . 

(63) 

If |u[’'^(to, 0 |l 2 becomes large or A > 0 is small enough compared to 
(to, ■)]^ 2 , then the second summand on right side of (I63|) is small 
compared to the first summand. Note that we may replace u[’'^(to,-) 
by multivariate spatial derivatives (to,.) for 0 < |/3| < m such 
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that an analogous estimate holds for spatial derivatives .)|^2 

for 0 < |/3| < m. We shall observe below that this small damping after 
discrete time (under the assumption that some exceeds 

a certain level (say 1) is strong enough in order to offset possible growth 
caused by the nonlinear terms. 

vii) We compare the damping with the upper bound of the nonlinear term. 

Remark 2.3. For analytical purposes the damping estimate ^ 

(for a short time interval length Ag) can even offset the growth of a rough 
upper bound ~ prAg (even the rougher upper bound ^ pcAg) as an upper 
bound for the possible growth cause by the nonlinear term as the factor 
r is arbitrary large such that pr^ can be chosen to be large compared 
to pr (especially such that for any given Ag pr^Ag is large compared to 
prAg or even compared to pr). However, we use the finer estimates above, 
where this may be also of computational interest. The offset of possible 
growth caused by the nonlinear term becomes effectives as the time interval 
Ag > 0 is small enough and the spatial parameter r > 1 (indicating 
deviation of the operator class from a strong semigroup property) is large 
enough. 

Summing up the preceding argument recall that we may replace r by r —tg 
on order to have for given to > 0 and t £ [to, + ^o] 

{t, x)\ < |i:)f'y[’‘'(to,-)L2exp(-47r2i/pr2 (T-tg)A2 ) 

+c^ [SDn'^pr'^jy^T - tg)Ag+^) + prLmi^npr'^v)^ ((r - to)~^~^C* + e) . 

(64) 

Recall from dMl) and analogous estimates for spatial derivatives that we 
have a damping estimate 

|^f<’'’’''(^o,-) *sp G^(r, .)|^2 

< (to, .)L. exp (-dTr^z^pr-VAg) + {8D^^pr^i.rAl+^) . 

(65) 

These estimates become effective for small Ag > 0 if the norm of (some 
spatial derivative) of the initial data exceeds a certain threshold, i.e., if 
for some (3 we have {to, .)|^2 ^ 1- If this threshold is realised for 

some 0 < |/3| < m, then the last upper bound term in (1641) or in (1651) is 
relatively small for small Ag > (t —tg) . Note that this upper bound term 
contains no spatial variables, i.e., it behaves like a constant with respect 
to spatial norms. Define 

Cd ■= Cn {SDir'^pr^iyTAl'^^) (66) 

For T £ [tg, tg + Ag] we have 

\D^'>J?’'"’''ito + Ao,.)\^ 2 < \vl'''ito,.)\^2exp{-ATT'^iypr'^){T-to)Al) + 
+prLjn{^'K'^pr‘^vY [Al~^C* + e) . 

(67) 
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It is sufficient to consider the case to = 0 as the same following estimates 
hold for to > 0 if T is replaced by r — tp. In this case the relation on (1571) 
shows us that 


( 68 ) 

if 

(exp(-47r2i/pr2(T-to)A§) - l) + 

(69) 

+prLm{^TT^pr‘^vY (Aj-^) C* + e < 0. 

Note that the positive real number e in (1531) is of exponential decay on a 
small time interval. More precisely, 

e ~ X)f=i /rd\b (ffi y)\d^dy i 0 as A ). 0, (70) 

hence e is comparatively small as A'Kpr‘^v{T — to) < dTrpr^i^Ao becomes 
small. Furthermore, as 0 as Ap ), 0 with A^+^, and the damping 
factor .)|^2 (l ~ exp (—dTT^z/pr^rAo)) is dominant for r = Ap 

if .)|^2 > 1 such that is relatively small compared to the 

modulus of the main part of this damping term 

|<’’'’‘'(0,.)L.(-4’rWAoA2). (71) 

Next we consider the conditions such that the modulus of the main damp¬ 
ing part is larger than the last term (1571) (the last term with factor e is 
comparatively small and can be neglected). Here we observe the expo¬ 
nents of the parameters p, r, v and Ap in dZH) compared to the exponents 
of the parmeters p, r and Ap of the last term in (1571) . For r = Ap for the 
damping term in d7T|) we have the dependence 

~ vpr^^l, (72) 

which has to be compared with the last term in (1671) . where we have the 
dependence 

(p)i+^(^)i+2^(^)^(Aj-^). (73) 

The estimate of convolutions of a Lipschitz function with first order spatial 
derivatives of the Gaussian forced us to assume lypr'^ > 1 which implies 
large r > 1 in general. Hence we have to impose 

5 i such that < r^. (74) 

Choosing a small step size parameter p, say p = Aq , we have 

(p)i+^(r)i+2'5Aj-‘5 = (75) 

which has to be compared with the damping term, where for p = Aq the 
latter has the dependence 

pr^Ag = A^+^r^. (76) 
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For 6 < 0.5 and small Aq (especially Aq < 1) we get the time step size 
condition 

2 “ 1 “ S 

fi{l + 5) + l — (5>^ + 3iff/r> —-— (77) 

0 

which implies ^ > 5 for 6 < 0.5. Hence any choice with r as above, 
(5 € (0, 0.5), and p as above with p > implies a regular upper bound 
for v"^,! < i < D. Note that r > 1 is related to a deviation from a strong 
semigroup contraction principle of the operator, and may be chosen large 
such that lypr^ > 1 holds for given z/ > 0. Note that in any case the upper 
bound constructed is at discrete times lAo,l € N, where N denotes the set 
of natural numbers. However, using the local contraction result and the 
semigroup property we get a regular upper bound for all time. 


3 Appendix A): Proof of local time contraction 


Again remember the notation discussed in assumption c) above. We consider 
the local contraction Lemma O in more detail. It is sufficient to prove the 
theorem for the functions v^{t,x) = where Vi^t{t,x) = v1^{t^x)j for 

pr = t. We then have 


= v,{to ,.) GP +pEf=i 

= v,{to ,.) GP +pEf=i Fl^{vP’>^-\VvP’>^-^) * GP, 


(78) 


where Gp = GP’^ in the notation above, i.e., r = 1. For 0 < I 7 I < m — 1 and 
related 1 < |/3| < m, where j3p = jp + 1 and Pi = 7 ; for Z 7 ^ p we have 


= D^,(to,.) GP + pEf=i D2FI^{vP’>^-\VvP’>^-^) * GP^p 
= DPvPh, .) *,p GP + pDlf\vP^’^-\yvP^>^-^) * Ge,p. 


(79) 


According to our notation discussed in item c) above, we have Fjj = fi for 
all 1 < i < D, so for fixed time r G [to,Zo + Aq] from ((751) we get 

= Eo<|/3|<m |-Df<’''(T,.) - .)|i2nc 


< p 
+p 


{fi{vP’’^-\VvP’^-^) - fi{vP’'^-\\7vP’^-^)) * GP 
{MvP’’‘-\^vP’’^-'^) - f,{vP’’^-'^,VvP’^-^)) * GP 


L^nc 


L^nc 


+ So<|7|<m-l P 
^ 0 <| 7 |<m —1 P 


{D2fi{vP'^-^,VvP'^-^) - D2f,{vP'^-^,VvP''^-'^)) * G0_^ 
{D2fi{vP’^-^,VvP’'^-^) - DlfPvP'^-‘^,VvP^^-‘^)) * GP^^ 


L^nc 


L^nc 

(80) 


If Fj has a local interpretation, then for all 0 < |/3| < m D^Fj is assumed to be 
Lipschitz continuous with respect to all arguments on compact domains, hence 
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DJf is also Lipschitz as a sum of Lipschitz functions. We may use an uniform 
Lipschitz constant L which serves for all arguments and get 


I P^k 

|< 

< p 
+p 






L\v', 


L\v‘. 


p.fc-i/ 




L^nc 


L^nc 


(81) 


+ Eo<|.|<™-i P .) - .)\* GP 


.fc-2/ 


L^nc 


+ Eo<|7l<»^-i P .) - .)| * GP 


.fc -2 / 


L^nc 


For a global operator Fj with Lipschitz continuity with respect to l-liSpn;;; norm 
a similar relation as in (IRT|l holds. We have considered this elsewhere in the case 
of the Navier Stokes operator, and may reconsider this in Appendix B below. 
As an essential case we estimate a summand of the third term on the right side 
of dSH). We write 


G'O.j — Isf (O)G'Oj + l(RD\Bf (O))G'0j, 


(82) 


where lilf’(O) is the characteristic function of the ball of radius 1 around 0. 
It has the upper bound 


L\D2vlf-\T ,.) - .)| * iBf (o)G^., 


L^nc 


+P 


.) - (r, .)| * ( 0 ))G^,, 


L^nc 


< pLC,^\Dlvt’}-\T ,.) - .) 




where 


G17 J^ \ lgpx)\dxdt + J ^l(RDy3D(Q))G0 , 


L^nc 


The convolutions with Gp can be treated similarly. We may define 
G07 •“ IIsf (o)GO(t, x)+ T ^l(RDy 5 D(Q))GP^ 


L^nc 


Hence 


”f'=(T, (r,.) 


p,fc-i. 




L^nc 


(83) 


(84) 


(85) 


( 86 ) 


+ Eo<|/3|<mP-^(Go7 +G 17 ) ^ {t, - DPv^’^ ^(t, .) 

We may then choose p = —rr^ - ^- 77 ; -;;—rr and get the desired contraction. 
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4 Appendix B): Additional remarks concerning 
local time contraction in case of the Navier 
Stokes equation 

We consider the local time contraction result in the special case of a global 
operator, i.e., in the case of the incompressible Navier Stokes equation in its 
Leray projection form. Here we may use that Lipschitz continuity of the Leray 
projection term function 


X —>• 



(87) 


holds, if 5 z, 1 < Z < D is located in a suitable strong function space, i.e., if 
gi € n C"" for m > 2 for all 1 < I < D. For coordinates (t, y) with y = rx 

and pr = t the Leray projection form of the Navier Stokes equation is 


dt 


D 

+ pr^ 
j=i 


dyP’^ 

dxj 



Vi'vz, 


{y)dy 


( 88 ) 

Remember: if we use viscosity damping in order to get global regular upper 
bounds, then we choose r > 1. This parameter can be used to measure the 
deviation from a strong semigroup property with respect to some strong func¬ 
tion space. A small parameter p > 0 is a useful tool in order to prove local 
time contraction. Next we reformulate the local time contraction result in this 
context. Assume that we have constructed a global regular upper bound for 
the original equation up to time to > 0. Now data Vi(to,-), 1 < i < D with 
Vi{to,-) G 77™ n C™ for some m> 2 are given. We define a time local iteration 
scheme 1 < i < D, /c>0ina time interval [^ 07^0 -I- A], where to > 0 and 


,,PX,0( 


(t,.) ='!;f’’'(to ,-)7 t€[to,to + A] (mb. i;°(0,.) =/i*(.)), 


(89) 


and is a solution of the linearized equation 


= PP /rD KdA- - y) < 


p,r,fc—1 p,r,fc—1 


(90) 


V. 


m,l 


(.y)dy, 


where 1 < i < D. This is a variation of the scheme above, where the first member 
at fc = 0 was defined by a convolution with a Gaussian. For a contraction in 
a 77™ n C^-norm (with respect to spatial variables) this variation of a scheme 
is sufficient (the result may be easily reformulated for the former scheme above 
and is the adapted for even the stronger function spaces considered). For the 
parametrized incompressible Navier Stokes equation the local contraction result 
in m has the foollowing counterpart. 

Lemma 4.1. Let to ^ 0 assume that for some m > 2 and a finite constant 
Co > 0 we have 
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Then there exists a A > 0 dependent only on dimension and on the constant 
\v^’^{to, ■)Isuc/i that on the time interval + A] the functional in¬ 
crements ^ < j < D satisfy 


sup 

rG[tojio+'^] 


\SvP 


,r,fe+l 


(a-) 




sup 

TG[to,to+A] 


IK” 


'(A.) I 




(92) 


and 

sup \SvP’^'\t,.)\ <^. (93) 

TG[to.to+A] ^ 

In order to prove this lemma (as a footnote to Lemma 12.11 and its proof in 
Appendix A) we first consider classical solution representations of the time local 
iteration scheme. We have 


P^T.k 

vr = v, 


r (to,.)*.pG^.’^ + pr(EK 


— 1 


dxi 


* G^’^ 


+p^ {iwn KdA- - y) EK=i y)dy) * 


(94) 


where we recall that G')’’’ is the fundamental solution of — pr’^G^’^ = 0, 
considered on the time interval t G [to, + A]. For 0 < I7I < to — 1 and related 
1 A |/3| < nT-, where /3p = 7p + 1 and Pi = ji for I ^ p we have 


P.r.fc ^ jjS 


\to,.)*sp GP/ 


prDl (E 


D P,r,k-l dvP'^'^ 
j=l Pj dxj 


-) 




+prDi (4^ KD,^{■ - y) Efm=i y)dy) * GO:;. 


Note that for fc = 1 we have 

c P,r,l P,r,l P,r,Q 

OV, = V, - V, 


(95) 


= ur(to,.)*.pGPE-K(io,.) 


+pr (Ef=i <’’'(^0, ■)) * G^’^ (96) 

+PP (/rO KD,^{. - y) Efm=l <m(^0, (to, y)^!/) * Gg’A 

Similarly, for 0 < I7 I < to — 1 and related 1 < |/3| < to, where /3p = 7p + 1 and 
Pi = 7i ioT I ^ p we have 

sDPvr^ = DPvr^ - DPvr° 

= DPvt^{to,.)*sp Gg’’'-i4^r’° 

/ \ (97) 

+prD2 (Ef=i K(io, OKf (^0,.)) * G^:;, 

(/rd A:D,,(.-y)EK=i<™(io,2/)K/(*0’y)‘^2^) * Kp- 

As |-Dfu((^’'(<o, ■)|i2p|(;7 (recall definition of the main text) spitting intergals in 
dMD and (1971) in a local part and its complement and using Young inequalities 
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and standard estimates it is straightforward to prove that given A > 0 and 
r > 0 there is p > 0 such that 


sup .)| < 

TG[to,to+A] 2 


(98) 


Next we may abbreviate 

^ prD2fr = PrDl (Sy 


D p^r,k— 1 dv^’' 


dxj 




+prD2 (v Kc,,(, - „) 


(99) 


and get 


p,r,k 






= E 


0<|/3|<m 


\DPvr\T,.)-DPvr’’^-\T,.) 


L^nc 


< pr 




L^nC 


L^nc 


Eo<\j\<m-i pr\ 

( 100 ) 

We estimate the last term in the latter equation (the estimation of the other 
terms is similar). For each 0 < jql < to — 1 we have 


pr 


= pr 


D2 (Ef=i * GZ + D2 (Ef=i 


L^nc 

-1 

dx-j 


+‘^ppDl Kdp{. - y) 


* G^’’' 

^ u,p 


L 2 nc 




( 101 ) 


where - according to our previous notation- — Note that 

for some m > 2 we assumed 






<Go 


for some finite constant Gq > 0. Assuming for /c > 1 


sup v: 

rCfioi^o+A] 


p,r,k — l 


(a-)| 




< Gfc_i 


( 102 ) 

(103) 


from (110411 and (IIOIII we get 

sup^g[t„_t„+^] I(5uf(r,.) 




< prcDmCcCKCk-i ^(r,.) 


(104) 




Here com is determined by the number of terms on the right side of (|94l) plus 
the number of terms (of derivative expansions) of (1551) for all 0 < |^| < to. 
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This constant may be determined by elementary combinatorics for numerical 
purposes (which is not our main interest here). Furthermore we may use local 
integrability of the Gaussian and first order spatial derivatives of the Gaussian 
on open time intervals (by standard Gaussian upper bounds, cf. main text 
above) and choose 


Cg Li((o,A)xR°) Si=l I Li((0.A)xR^) 

"^Sp=l I lR-D\Bf ( 0 )^ 1 )’ Ilp((o,A)xR-°) |lR°\Bf’(0)^j/j'| lp((o,A)xR-°)' 

(105) 

As for the constant Ck we may also use local intergability of the first order 
derivatives of the Laplacian kernel and L^-integrability outside a ball and choose 

Ck = + Si=l 

(106) 

+ |lR°\Bf (0)^G',iL2((o_A)xR^)- 


Next note that 

+ EtlSUp^G[to.to+A] 


(107) 


such that the choice 


leads to 


1 

^ “ 2rCDmCGCK{Co + 1 ) 


(108) 


sup |<’’'’''(t, ■)\jjrr^nc^ < Go + 1, (109) 

re[to,to+A] 

or Ck-i < Go + 1 for all /c > 1 in the induction hypothesis in (I108|) above. 


5 Appendix C): Additional remarks concerning 
randomized models (appearance of supercrit¬ 
ical barriers) from the perspective of Trotter 
product formulas 

Finally, we remark that additional ideas are required if this method is to be 
applied to randomized models, e.g., random averaged Navier Stokes equations 
or Navier Stokes equations with stochastic white noise. The reason is that the 
viscosity damping can only offset possible growth caused by the nonlinear terms 
if the functional increments of the Burgers term and the Leray projection term 
(considered over a small time interval + A] with A > 0 arbitrary) live in 

strong spaces. For example, if the problem is posed on a torus (periodic bound¬ 
ary conditions), and regularity is measured by the order of decay of the modes 
of a local solution function as frequencies become large, then the decay of the 
functional increments of the Burgers and Leray terms is determined essentially 
by the decay multiplication of infinite matrices of modes applied to the data 
(with a diagonal weight for the first order derivative of the Laplace kernel in 
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the case of the Leray projection term). If a noise term enters the equation or a 
random averaged model is considered then a direct application of the iteration 
scheme would lead to the effect that modes are mixed up and even properties 
which may be preserved pathwise may not be preserved ’in expectation’. Depen¬ 
dent on the model of course, supercritical barriers appear for strong solutions, 
and stronger damping would be needed in order to establish global existence of 
strong solutions. Let us consider this in more detail and consider the Navier 
stokes equation problem posed on a torus of dimension D and diameter I 
(the diameter may be set to Z = 1 of course). 

Given a spatial scaling parameter value r > 0 and writing the velocity com¬ 
ponent vl = vl (t, y) for fixed t > 0 in the analytic basis {exp ^ q, g 2^} 

vl{t,y)-.= ^ Via{t) exp (hq) 


the spatially scaled equation in (1^ and (1^ becomes an infinite ODE equa¬ 
tion for the infinite time dependent vector function of velocity modes vl^, a € 
iP , 1 < i < D of the form 


dt ~ P^j=l ) ^ia ^ 2^j=l 2.^7eZO 

+rZ7rtailfa^0} ' 


27r27j 


( 111 ) 


for all 1 < f < D and where for all a £ iP we have f[„(0) = hia along with 
the a-modes hia of hi. In a calculus with infinitesimals, we may describe the 
solution by Trotter product formulas obtained form iterative infinitesimal Euler 
steps. As an example consider nonstandard calculus (we do not enter techni¬ 
cal peculiarities and describe ’internal sets’ and ’extension of a set-theoretical 
universe’ as these specifications are rather trivial or well-known at least). For 
arbitrary tg > 0, a hyperfinite number N and an infinitesimal 6t with N6t = te 
we may define a nonstandard scheme for time steps 


mSt £ {0, St, 2St, • • • , {N — l)St, NSt = 7} . 


We have the infinitesimal Euler scheme 

+ l)St) = vlaimSt) + ^ Via{mSt)St 




St, 


or with the abbreviation 


elja-yi'mdt) = -r 


27rj(o!j-7j) 




+r2mail{a^Qp'K' 


2 Efc=i 


{rn5i) 


I 


( 112 ) 


(113) 
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we have 


<a((w + l)St) = v^JmSt) + X;f=i vl^{m5t)5t 

(114) 

+ Ef=i elja-fim5t)v^j^imSt)dt. 

In the case Z = 1 we have a solution representation 

exp (-J^r 247 r 2 J^Zi x 

(115) 

Here the symbol = means that the identity holds up to an infinitesimal error, and 
we denote v’’’^ = (v[’^, • • • )'^ with D infinite vectors vl'^ = ('c[„)aGZ° • The 

entries in {Sijap) are Kronecker-Js which describe the unit DZ^ x UZ^-matrix. 
(We note that internal counterparts of the letter sets should be considered in 
order to draw global conclusions using the usual methods- again these extensions 
are rather trivial). The formula in (11151) is easily verified via formulas at each 
time step m of the form 


|^(5„a;3exp (exp (((eCj„^(m5t)) <5t)) v^’^imSt) 

(116) 

(as a representation for v^((m-|- l)(5t)). The use of explicit infinitesimals allows 
us to have an effective use of first order equality = for arbitrary finite time where 
on an infinitesimal time level we have simplifications of the formula in (11151) in 
the sense that 


V 


\te) = (%a/3 ( 1 - i^47rV^a2(5t j j (l -t- (mSt)^ 6tj h^, 


(117) 

is also valid up to order 0{St^) (if 6 t is infinitesimal). Similarly, at each time 
step number m an infinitesimal Eulerstep (with error 0(6t^)) is described by 


JlJOl) 


I exp 


-z/dTrV^ a^Jt \ ) (l + [{<jc.p{mSt))Stj v^-^{mSt). 


(118) 

We observe that at each time step all modes get an viscosity damping except 
for the zero modes. However this zero modes can be controlled and are a minor 
problem. We considered several ways do deal with this. Here we consider an 
extended Trotter product formula where the zero modes are shifted to zero at 
each time step and reappear in the drift of an equivalent related problem. This 
is effectively a control function for the zero modes (we use a different notation 
here than we did elsewhere). 

Define , ■ ■ ■ vW’^ Y with D infinite vectors u)" 


c,r,F / c,r\ 
= )c 


/ c,r,F c,r,F\T 

(^1 ! ■ ■ ■ '.'j — \'^ia Jai 

where at all time steps m > 0 we have (mSt) = vl^ (mSt) — Cia (mSt) where 
Cia = 0 for a ^ 0 and Cio{mSt) = —vlf^irnSt). We may plug in vl^(m5t) = 
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''^ia + Cia {mSt) = (mSt) into the equation in (I112D and get 

dt ~ 2^3 = 1^^ P j^ia ^Z^j=lZ^7GZ^ i \^^j(a-7) + ^{a=7}''j(a-7) 

E”fc=i E^gzO 47r^7j(afc-7fc)^^°;''<’(I,_^) 


-t-r27r2c^jl|Q,:^Q| - 


Ef=i 47r2c 


(119) 

for all 1 < 2 < D and where for all a € iP ^ and where we have {0) = hia 
for a 0 and 2;((^’^(0) = 0 for all 1 < 2 < D. In the case / = 1 we have a solution 
representation 


= n^Jo (dya/3 exp '£?=! <4^)) X 




X I exp 


( 120 ) 


where h'^’^ = 2 ;'^’’'’'^( 0 ) are obtained form the data by setting the zero modes 
to zero. All related formulas for uncontrolled schemes above have their obvious 
counterparts which are obtained by substitution of I’y Ob'^ioeisly 

we have 


+r27r2aa{„^o}47r^--, 


( 121 ) 


Ef=i4^"c 


Now assume that the data hi are in a strong Sobolev space, lets say hi £ 

fjD+2 

in case of dimension D. The dual Sobolev norm 


h,eh"^ = /i™ (Z^) iff ^ |/2,cP(l + lap"") < oo 


( 122 ) 


tells us that there surely is a finite constant such that 

VaeZ^: |/ 2 ,„| < <oo. 

v'=’’'’^((to + l)6t) = [hijap exp {-vAPP YJiLi x 

Recall our notation; v°’’'’^(m5t) = (mSt), ■ ■ ■ , ^^"^(mdt)), where up "^(rndt) = 

In the Euler scheme (one Euler step) the Euler terms can be estimated by 
elliptic integrals, such that for some constant c (depending on dimension D) we 
have 


(123) 

(124) 


27r{D + D^) J2 


me 


c 


1 + |a - /3P+2 I + |flP+2 ^ 1 


cC^ 


a 


D+3 ■ 


(125) 
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i.e. the order of polynomial decay is preserved (even more than preserved). 
Further transformation to another time scale t = pr may be used together with 
a strong damping parameter r > 1. The nonlinear terms then get the scaling pr 
while the viscosity coefficient (the coefficient of the Laplcacian) becomes pr^v. 
For given v > Q the parameter p may be chosen (along with large r > 1) such 
that this latter term pr'^v may be large while pr is small such that the scaled 
upper bound variation of (11251) becomes 


pr2'K{D + D"^) ^ 


|/3|C 


C 


< 


prcC"^ 


1 + |a -/3P+2 1 + I/3P+2 - 1 + 


|r>+3 ■ 


(126) 


where prcC^ become small such that geometric series upper bounds for the sum 
of functional increments can be constructed. We have discussed such niceties 
elsewhere, and we do not need them if we want to construct an upper bound for 
fixed T > 0. For the regularity of the data we assume that there is a constant 
C > 0 such that for all 1 < i and all a S lP we have 


hiry ^ 


c 


■ - l + |aP+2- 

Given any T > 0 and data as in dUZI) we may choose (generously) 

1 


r > - {C + cC^il+T)) 


(127) 


(128) 


and get a finite constant C* = (G + cG^(l + T)) such that for all 1 < i < H 
we have a global regular upper bound 


te [0,T]. (129) 


This theory cannot be applied to stochastic versions of the Navier Stokes equa¬ 
tion, however- due to the fact that elliptic integral relations as in (112511 are 
sensitive to such extensions. As an example consider a Navier Stokes equation 
model with a stochastic force (white noise), i.e., an initial value problem of the 
form 


dt 


D 

vAvi + ^ 
i=i 



„ dW 


(130) 


along with incompressibility J^iLi = 0 and data Ui(0, .) = hi G fl G'" 
for some m > 2. Here the derivative of the Wiener process W may be defined 
via the represention 


W{t) = YyNn 

n—1 


s)ds, 


(131) 


where (lV„))?Li is a family of independent, identical, standard Gaussian ran¬ 
dom variables, and denotes a family of orthogonal functions define for 

nonnegative time (this orthogonal family is orthonormal in the case of standard 
white noise). We have 

oo 

— = Y,Nnkn{t). (132) 
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An extension of the Trotter formula scheme above (on an infinitesimal time 
scale) leads to an additional source term 

OO 

W{mSt) — W{{m — l)St) = ^ Nn{kn{mSt) — — l)St)). (133) 

n—1 

For such a scheme energy can be transported from lower to higher frequencies, 
and -depending on the family (A:„) upper bound relations of the form (11251) 
may not hold for the solution. Similar consideration apply to stochastic average 
models. 


6 Appendix D): Extensions of the result to mod¬ 
els with generalised diffusion term 


The Gaussian damping estimate in item vi) can be generalised to models with 
variable diffusion. Such generalisations are desirable as even hydrodynamic 
limits of the Boltzmann equation around the Maxwell state lead to a viscosity 


ly 




(134) 


Here D* is the Legendre dual of the Dirichlet form of the collision term of the 
Boltzmann equation linearized at the Maxwell equilibrium. For models with 
uniform elliptic diffusion terms, i.e., terms of the form 


D 


E 
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dxjdxk ’ 


where for some constants 0 < A < A < 00 


(135) 


D 

y{t,x) G [0,oo) X ^ < A|^p, (136) 

i,fc=i 

there are Gaussian upper bounds of the fundamental solution. The estimates 
of item vi) can the be adapted to this situation straightforwardly using these 
Gaussian upper bounds. 

For generalisations to highly degenerate diffusions we have to recall Horman- 
der’s result. Let us For positive natural numbers m,n a matrix-valued function 


X —>• {vji)'^'^{x), I < j <n, 0 < i <m, (137) 

on R", defines m -|- 1 smooth vector fields 


i=i 


d 

dx/ 


(138) 


where 0 < * < m. Hormander provided a sufficient condition for the fundamen¬ 
tal solution on [0,oo) x R" of 


^ = ^EZiK^p + VoP 
p{0,x;y) = Sy{x), 


(139) 
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where Sy(x) = S{x — y) is the Dirac delta distribution with an argument shifted 
by the vector y € R". The condition is: for all x G K” assume that 

H, = R", (140) 


where 

:= span|D,(a;), [Vj-, 14] {x), [[V, , 14] , Vi] {x), 

(141) 

•••|l<i<rn, 0 < k,l, ■ ■ ■ < rn^, 

and where .] denotes the Lie bracket of vector helds as usual. Assume that the 
coefficients of the vector fields are smooth (i.e. C°°) and bounded with bounded 
derivatives, i.e. Vji € (7“ (R”) (sometimes linear growth for the functions vji 
themselves is allowed ( Vji G C^i (R") in symbols). Indeed this makes no real 
difference. These strong regularity assumptions are due to the fact that i Then 
a smooth density p exists and for each nonnegative natural number j, and 
multiindices a, /3 there are increasing functions of time 

Bj,a,p • [0)2^] (142) 


and functions 
such that 


'nj,a,p, : N X X —>• N, 




(143) 


(144) 


< 




exp (-5j,a,/3(t)^^7^) • 


As the damping estimate is in discrete time at each time step of the scheme 
above the local damping estimate above in item vi) can be adapted to this 
situation. This is a first step in order to extend the construction scheme of 
global regular solutions models of the form 


^ 1 Eti V?v, + Vov, + /.(u, Vu) = 0, 1 < * < D 

Vi{0, .) = h^. 


(145) 


A stronger form of condition d) is required in order to get an upper bound for 
the functional increments 

n pto-\-A p 

6vi{t,.) =Vi{t,.)- v^{to,.)p{t,x]s,y)dy = j / fi{v{s,y),Vv{s,y))p{t,x;s,y)dyds, 
JS.P' Jto Js.P’ 

(146) 

because we have to compensate for the factor (l+a;)"*^'“''^ in the Hormander esti¬ 
mate of the density p (which satisfies (11451) 1. Note that some order of spatial de¬ 
cay is lost in general as the (smoothed) initial data term Vi{to, ■)p(t, x; s, y)dy 
shows. However this loss of spatial decay is not increased iteratively if the non¬ 
linear term satisfies a strong submultiplicative property which compensates the 
factor (1 -Px)™^’'*''^. Clearly, the loss of polynomial decay depends on the order 
of regularity in terms of the constants rrij, but such a loss is not iteratively 
increased by the linear term as a Chapman Kolmogorov semi-group property 
holds for the density p and since the functional increments preserve polynomial 
decay for an appropriate strong submultiplicative property. 
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7 Conclusion (and comparison with results form 
CKN theory) 

The preceding argument shows that there is a considerable class of partial dif¬ 
ferential equations (subsuming the incompressible Navier Stokes equation) with 
spatially subhomogeneous nonlinear terms which have global solution branches 
in strong spaces for strong data. This conclusion contradicts the view of many 
authors that global regular solution branches of the Navier Stokes equation 
have to be due to the special structure of the equation (especially the Leray 
project term) if there are any such solutions. A difficulty of the construction 
of global regular solution branches for some members of some classes of equa¬ 
tions arises if the equations are considered in weak function spaces, i.e., for data 
which live in weak functions spaces, for equations with time-dependent external 
forces, or for equations with stochastic force terms or stochastic average equa¬ 
tions. Especially, in the latter cases the strong upper bounds of weakly elliptic 
intergals or related multiplication rules of infinite mode matrices do not hold 
anymore for iteration schemes in general. Energy can then be transported from 
lower frequencies to higher frequencies and supercritical barriers may appear, 
i.e. models with stronger viscosity damping terms may be necessary in order to 
get global regular solution branches or, at least, to avoid the appearance of sin¬ 
gularities. It is well known that in the case of the incompressible Navier Stokes 
equation global regular solution branches are unique (this is a difference to the 
Euler equation as we argued elsewhere). If a global regular solution branch 
Vi € C° ([0,T],i7"‘ n C™), 1 < * < Z? for arbitrary T > 0 is given in the case 
of the Navier Stokes equation, then a Cornwall inequality implies uniqueness, 
i.e., if Vi, 1 < i < D is another solution of the incompressible Navier Stokes 
equation, then 

\vit) - v{t)\l^ < \v{0)-v{0)\l^exp(^C (lh('S)lL4 + ||i^(s)|i4) (147) 

for some p > 4 which depends on the dimension (here p = 8 in case of dimension 
U = 3 is sufficient), and where C > 0 is a constant which depends on the 
dimension D and the viscosity. Here the function space fl C™ is appropriate 
for Cauchy problems on the whole domain, where these function spaces reduce 
in the case of a torus, of course. Such a Cornwall inequality implies that a 
solution branch is unique in sufficiently strong function spaces. However we 
have no such inequality for the whole class considered here. The existence of 
global regular solution branches depends on a local time contraction theorem in 
strong spaces. Such local time contraction results in strong spaces can also be 
used in order to obtain global regularity and existence results as a consequence 
of the CKN result. In this case it is essential to show that a Leray-Hopf solution 
is left continuous with respect to time and in strong function spaces with respect 
to the spatial variables. If tp x is a time slice where a potential singularity 
of the Haussdorff set of singularities is located, then we have argued elsewhere 
that a Leray Hopf solution Vi, f <i < D, which is smooth on [tp — e, tp) is left 
continuous at tp such that 


Ui(tp,.) G n C''" ,f < i < D, for some m>2. (148) 
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This is the essential step of course as we have local time contraction and can 
combine the standard techniques of construction of a Leray Hopf solution with 
local time contraction and the singularity analysis. Having obtained (11481) we 
may obtain a local regular solution on a time interval [to, to + A] by local time 
contraction (cf. also Theorem 17.II for an alternative local time result) (we may 
then start with this local solution and extend it applying the standard techniques 
of projection smoothing and the usual limit procedures. The resulting extended 
solution has a thins set of possible singularities at some time sections ti x for 
ti G I in an index set I of Lebegues measure zero and ti > to+A. For the smallest 
such ti we may consider a singularity analysis -as we have shown elsewehere- 
and get a regular solution on a time interval [to,ti] with Vi{ti,.) € D C™ 
which is left continuous at t/. We may then apply the local contraction result 
for the new data .by the semigroup property- and repeat the procedure). Here, 
we can use the local time contraction results in n C"* norms (spatially). 
Note however, that we can also use local time regularity results such as 

Theorem 7.1. For some to > 0 assume that Vi{to, ■) G ,1 < i < D, p > D, 
^^^iVi^i{to, ■) = 0. then there is A > 0 and a unique short time solution 
i>i, 1 < I 'A D on the time interval [to, to + A] such that 

v^gC ([to, to + A], LP (T^)) n ((to, to + A) X T^) (149) 

for all 1 < i < D. 

Hence it seems that global regularity results can be obtained for data (T^) 
and p > D ii the problem is posed on the torus. Note that this also implies 
that regularity implies global smoothness. This is due to the well-known 
embedding 

HP,s ^ foj. g > p, s > and i i 4 (s - t). (150) 

p q D 

Here is the usual analog of For example, in case D = 3, p = 2 
and q = A we get n 774 , 0.25 ^ ^4 Similarly for H = 3 and g = 3 -I- e for 
small epsilon we get an embedding in Hence it seems that even data 

are sufficient. This results certainly is stronger then the results of this paper, 
but the CKN results cannot be transferred to the large class of equations which 
is considered here. This class (and variations of it with a more general class of 
viscosity terms) may be important also in a physical sense, as data imply that we 
have to go beyond the simple Navier Stokes equation in order to describe fluids 
appropriately. Furthermore the approach considered here is constructive and 
the regular upper bounds constructed can be a guide for algorithmic solution 
schemes. Finally we note that the existence of a (possible highly degenerate) 
diffusion term is essential for the preceding argument. Viscosity limit arguments 
[v 4 , 0 ) require additional specific structure (such as vorticity forms in the case 
of Navier Stokes and Euler equations), and even then limits are of restricted 
regularity in general. 
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